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Abstract 

> ' 

Elaborating on a previous work by Han et al., we give a general, basis- 
independent proof of the necessity of negative probability measures in order 
' for a class of local hidden-variable (LHV) models to violate the Bell-CHSH 

inequality. Moreover, we obtain general solutions for LHV-induced probability 
measures that reproduce any consistent set of probabilities. 
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In 1982, Muckenheim [1] made use of negative probability functions in an at- 
tempt to resolve the Einstein- Podolsky- Rosen (EPR) paradox [2]. For this purpose, 
Muckenheim built a classical model endowed with negative probabilities that re- 
produced all the statistical predictions of quantum theory for the singlet state of 
two spin-half particles. While the physical meaning of extended probabilities is 
far from obvious [3], this attempted solution of the EPR paradox might seem [4], 
". . . as unattractive as {but not more unattractive than) all the others" (present au- 
thor's emphasis). Subsequently, Home, Lepore and Selleri [5] put forward a general 
argument demonstrating that one can always reproduce the quantum mechanical 
results for nonfactorizable state vectors of correlated systems by means of probabil- 
ities of the Clauser-Horne type [6] provided one allows for probabilities not obeying 
Kolmogorov's axiom according to which probabilities p are restricted to the range 
< p < 1. More recently, Han, Hwang and Koh [7] obtained explicit solutions for 
classical probability measures that reproduce quantum mechanical predictions for 
some spin-measurement directions for all entangled states, and proved the necessity 
of negative probability measures in this case. In the present paper, we shall extend 
the proof by Han et al. in the following sense. While Han et al.'s proof relies in 
a special basis (see Eqs. (28) and (39) of [7]) to show the negativity of probability 
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measures for the considered local hidden- variable (LHV) model, ours proves the ne- 
cessity of such negative probability measures in all instances where the predictions 
of the LHV model are made to violate the Bell-CHSH inequality [8-10]. This is done 
without relying on any particular basis states or measurement directions. In fact, 
our result holds irrespective of whatever quantum mechanical consideration. More- 
over, we give general solutions for LHV probability measures that reproduce any 
conceivable set of probabilities satisfying certain requirement conditions, namely, 
the normalization condition and the causal communication constraint (cf. Eqs. (17) 
and (28)- (29) below). A set of probabilities fulfilling these requirement conditions 
will be referred to as a consistent set. The proof goes in a rather straightforward 
way as follows. 

Consider an experiment of the EPR type designed to test the Bell-CHSH in- 
equality. Two correlated particles 1 and 2 fly apart in opposite directions from some 
common source. Subsequently, each of the particles enters its own measuring ap- 
paratus which can measure either one of two physical variables at a time — a>\ or a 2 
for particle 1 and b\ or b 2 for particle 2. The possible values of these variables may 
be taken to be +1 and —1. The source emits a very large number of particle pairs. 
The basic entity to be considered is the joint probability p(aj = m,bk = n) that the 
outcome of the measurement of aj on particle 1 is m, and that the outcome of the 
measurement of bk on the paired particle 2 is n, where j, k = 1, 2, and m,n = ±1. 
A representative deterministic LHV model describing this experiment is as follows 
[7,11]. The main assumption made by such a model is that, for every pair of particles 
emitted by the source, there exists a hidden variable A (with domain of variation A) 
which determines locally (for example, at the common source) the response of the 
particles to each of the measurements they can be subjected to. For the experiment 
under consideration, the set of all A can then be partitioned into 16 disjoint subsets 
Aj (with respective probability measure m«) according to the outcomes of the four 
possible measurements, a± and a 2 for particle 1 and b\ and b 2 for particle 2. In 
Table 1 we display the 16 rows characterizing the subsets Aj. The ith row indicates 
the response of the particles to the different measurements when the particle pair 
is described by a hidden variable pertaining to the subset Aj. So, for example, if a 
particle pair is described by a given A G A 2 , then the particles must behave accord- 
ing to the following local plan: if a\ is measured on particle 1 the result will be +1, 
if a 2 is measured on particle 1 the result will be +1, if 61 is measured on particle 
2 the result will be +1, and if 62 is measured on particle 2 the result will be — 1. 
(Note that, for each of the plans, the agreed result for aj is independent of which 
measurement (b\ or 62) is performed on particle 2, and similarly the agreed result 
for bk is independent of which measurement (ai or a 2 ) is performed on particle 1.) 

From Table 1, we can readily compute the predictions that our LHV model makes 



2 



Table 1: The 16 possible subsets into which the total set A can be partitioned. The hidden 
variables in each subset Aj uniquely determine the outcomes for each of the four possible single 
measurements a\, bi, 0,2, and 62- 



Subset of A 






a 2 


b 2 


Probability measure 


Ai 


+ 


+ 


+ 


+ 


mi 


A 2 


+ 


+ 


+ 


— 


m 2 


A 3 


+ 


+ 


— 


+ 


m 3 


A 4 


+ 


+ 






m 4 


A 5 


+ 




+ 


+ 


m 5 


A 6 


+ 




+ 




m 6 


A 7 


+ 






+ 


m 7 


A 8 


+ 








m 8 


A 9 




+ 


+ 


+ 


m 9 


A10 




+ 


+ 




m 10 


An 




+ 




+ 


mn 


A12 




+ 








A13 






+ 


+ 


m 13 


A14 






+ 






A15 








+ 


77115 


Ale 










m 16 



for the various probabilities p(aj = m,bk = n). These are given by 



Pi 




p(ai+ 


bi+) 


= mi + m 2 + m 3 + m 4 , 


P2 




p(a>i+ 


bi-) 


= m 5 + m 6 + m 7 + m 8 , 


P3 




p(ai- 


61+) 


= m 9 + mio + mu + mu, 


Pa 




pifli- 


bi-) 


= m 13 + m i4 + mi 5 + mi 6 


Pb 




p(ai+ 


b 2 +) 


= mi + m 3 + m 5 + m 7 , 


Pe 




p{ai+ 


b 2 -) 


= m 2 + m 4 + m 6 + m 8 , 


Pi 




pifli- 


b 2 +) 


= m 9 + mn +m i3 +m 15 , 


Ps 




p(ai- 


62-) 


= m 10 + m 12 + m u + mi 6 


P9 




p(a 2 + 


bi+) 


= mi + m 2 + m 9 + mi , 


Pw 




p(a 2 + 


bi-) 


= m 5 + m 6 + 77113 + mu. 


Pll 




P(a 2 - 


61+) 


= m 3 + m 4 + mn + rn 12: 


P12 




p{a 2 - 


61-) 


= m 7 + m 8 + m 15 + m 16 , 


Pl3 




p(a 2 + 


b 2 +) 


= mi + m 5 + m 9 + mi 3 , 


PlA 




p(a 2 + 


b 2 -) 


= m 2 + m 6 + mio + m u , 


Pl5 




p{a 2 - 


b 2 +) 


= m 3 + m 7 + mn + rn 15: 


Pw 




P(a 2 - 


62-) 


= m 4 + m 8 + m i2 + mi 6 , 



(1) 

(2) 
(3) 
(4) 
(5) 
(6) 
(7) 
(8) 
(9) 
(10) 

(11) 
(12) 

(13) 

(14) 

(15) 

(16) 



in obvious notation. We are assuming throughout this paper ideal behavior of the 
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measuring apparata and, in particular, perfect efficiency of the detection equip- 
ment. This requires the probabilities p(a,j = m,bk = n) to satisfy the normalization 
condition 

^2 P( a j = m A = n) = 1, (17) 

m,n=±l 

for any j, k = 1,2. From Eqs. (l)-(4), this in turn implies 

16 

£>i = l. (18) 

i=i 

For convenience for what follows we define the following two quantities, 

Si = m 4 + m 5 + m 6 + m 8 + m 9 + mn + m i2 + mi 3 , (19) 

and 

S 2 = 1 - Si = m l + m 2 + m 3 + m 7 + mio + m M + m 15 + rai 6 . (20) 

Let us now consider the sum of correlations 

A = c(ai, 6i) + c(ai, 6 2 ) + c(a 2 , &i) - c(a 2 , 62), (21) 

entering into the Bell-CHSH inequality [8-10], |A| < 2, with the correlation coeffi- 
cient c(a,j,bk) being given by 

c{aj,b k ) = p(aj = l,b k = 1) +p(% = -1,6*, = -1) 

- p( % = 1, b k = -1) - p( aj = -1, 6 fc = 1) . (22) 

Substituting this in Eq. (21), and taking into account the normalization condition 
in Eq. (17), the quantity A can equivalently be written in the form 

A = 2(pi + p 4 + p 5 + p s + p 9 + P12 + Pu + P15 ~ 2), (23) 

where we have used the abbreviated notation introduced in Eqs. (1)-(16). Now, 
by replacing the probabilities appearing in Eq. (23) by their respective expressions 
in Eqs. (1)-(16), one obtains the prediction that our LHV model makes for the 
Bell-CHSH sum of correlations, 

A LHV = 2(1-2S 1 ) = 2(2S 2 -1). (24) 

The Bell-CHSH inequality is violated whenever |A| > 2. Then, from Eq. (24), it 
follows at once that in order for the LHV model to give a violation of the Bell-CHSH 
inequality it is necessary that either S x < or S x > 1 (or, correspondingly, that 
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either S2 > 1 or S2 < 0). In any case, to have |Alhv| > 2, it is necessary that 
either Si or S 2 be negative. Thus, as the negativity of either Si or S 2 implies the 
negativity of at least one of the m/s, we have proved the necessity of negative prob- 
ability measures for the considered LHV model if this is to violate the Bell-CHSH 
inequality. On the other hand, from Eq. (24), it is also clear that the prediction by 
the LHV model does satisfy the Bell-CHSH inequality if, and only if, the condition 
< Si,S 2 < 1 is satisfied. Needless to say, the usual case where J2 m i = 1 an d 
< rrij < 1 fulfils this latter condition, and hence, for such a case, |A LH v| < 2. This 
shows that the role played by the assumption of nonnegativity of the probability 
measures in deriving the Bell-CHSH inequality is as fundamental as the role played 
by the assumption of local realism. It should be noticed, however, that, although 
necessary, the requirement of negative probability measure is not a sufficient condi- 
tion in order to have |A L hv| > 2. Indeed, it may be the case that some of the mj's 
be negative, and yet having that < Si, S 2 < 1. 

Now we are going to give general solutions for LHV probability measures that 
reproduce any consistent set of probabilities {p±,p 2 ,... ,pw}- To this end, let us 
first consider the special case where A = 2\/2. This value for A can be reproduced 
by the LHV model if S : = (1 - v / 2)/2 and S 2 = (1 + y/2)/2. An immediate, 
"equally-distributed" solution of these equations is 

1 — \/2 

m 4 = m 5 = m 6 = m 8 = m 9 = mu = m 12 = m 13 = — — — , (25) 

and 

mi = m 2 = m 3 = m 7 = m w = m 14 = m 15 = m 16 = — — — . (26) 

16 

Substituting these values in Eqs. (1)-(16) yields the following positive probabilities 
predicted by the LHV model for the special case considered, pi = P4 = p§ = p$ = 
P9 = P12 = Pu = P15 = (2 + V2)/8 and p 2 = p 3 = p 6 = p 7 = Pw = p u = p 13 = Pl6 = 
(2 — \/2)/8. It is to be noted that these values for p 1: p 2 , . . . ,p\e are the same as those 
predicted by quantum mechanics (QM) in the case that the quantum Bell-CHSH sum 
of correlations attains the Cirel'son limit Aqm = 2\/2 [12]. It should be added that, 
however, the equally-distributed solution, namely that for which m 4 = m 5 = m 6 = 
m 8 = m 9 = m u = m 12 = m 13 and m\ = m 2 = m 3 = m 7 = m 10 = m 14 = m 15 = m 16 , 
is clearly too restrictive since, as may readily be checked from Eqs. (1)-(16), it 
invariably leads to the prediction that p\ = p 4 = p$ = p 8 = p 9 = p 12 = p u = p 15 
and p 2 = p 3 = pq = p 7 = pw = pn = p± 3 = pi$. Consequently, except for the case 
where these conditions on the probabilities are met, the equally-distributed solution 
cannot account for the generic set of probabilities {pi,p 2 , ■ ■ ■ ,Pie}- Anyway, it is 
nevertheless important to realize that, if such probabilities Pi,p 2 , ■ ■ ■ ,Pie are to be 
given by the LHV predictions on the right-hand side of Eqs. (1)-(16), then the 
probabilities themselves must obey certain requirement conditions. Specifically, if 
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we want the generic probabilities pi,P2, ■ ■ ■ ,Pie to be cast into the form displayed by 
Eqs. (1)-(16), then the following relationships between them must necessarily hold 



P2 



[1 - pi - p 4 + p 5 - p 8 - p 9 + p 12 + P14 - P15, 



p 3 = -(1 - pi - p 4 - p 5 + p 8 + P9 - Pl2 - PU + Pl5) 



1 +Pl -Pi -p 5 



P9 +P12 + Pu -P15) 



P6 = 

P7 = ^(1 - Pi + Pa - P5 - Ps + P9 - P12 - Pu + P15) 



2 

1 . 

PlO = ^-Pl+Pi+Pb-Ps-P9-Pl2+ PU ~ Pl5) 
Pll = ^(1+Pl-P4-P5+P8-P9-Pl2~ PlA + P15) 



Pl3 = ^-Pl+P4+Pb-P^+P9-Pl2- PlA ~ P15) 
1 . 

Pl6 = ol 1 +Pl _ P4 ~P5 +P8 -P9+P12 -Pl4 -P15) 



(27) 



We note, incidentally, that the set of conditions in Eq. (27) is equivalent to the 
conjunction of the normalization condition in Eq. (17) and the so-called causal 
communication constraint [13]. This latter consistency condition requires that [14] 

piaj = m,bi = n) = p(aj = m, b 2 = n), (28) 

ra=±l ?i=±l 

^ p(ai = m,b k = n) = ^ p(a 2 = m,b k = n), (29) 



m=±l 



m=±l 



for any j, k =1,2 and m,n = ±1, and prevents the acausal exchange of classical 
information between the two parties involved in the EPR experiment. Both quantum 
mechanics and our LHV model satisfy the requirement of causal communication and 
hence the predictions by such theories do satisfy each of the constraints in Eq. (27). 

In searching for a general solution for LHV probability measures mi,m2, . . . , 
m 16 that reproduces the probabilities Pi,P2, ■ ■ ■ ,Pw, what really matters is the fact 
that, as can be seen from Eq. (27), only eight of such probabilities are independent. 
We may take the independent probabilities to be pi, p 4 , p 5 , p 8 , pg, P12, Pu, and 
P15. Therefore, to invert the set of Eqs. (1)-(16), it suffices to consider the eight 
equations (1), (4), (5), (8), (9), (12), (14), and (15), plus the normalization condition 
in Eq. (18). We are thus left with a system of 9 equations with 16 unknowns 
mi, m2, • • • , in which the probabilities p±, p 4 , p$, p$, pg, pu, pu, P15 are treated 
as given parameters. This system determines 9 probability measures as a function 
of the remaining 7 probability measures and the 8 probabilities above. So, for 
example, we can get the following general solution for which the set of probability 
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measures {mi, 777.4, 7776, m 8, m 9, mn, 777.12, m.13} is given in terms of the remaining 
set {777 2 , 7773, 7777, 777io, m.14, mis, mi 6 } an d the eight probabilities p±, p 4 , p 5 , p$, p 9 , p± 2 , 



Pu, and pi 5 , 

777! = -( — 1 — 2777 2 — 2777 3 — 2777 7 — 2777i — 2777 14 — 2777i 5 — 2777i 6 

+ Pi + Pa + P5 + + P9 + P12 + Pu + P15), (30) 

1 

777 4 = -(1 + 2777 7 + 2777 10 + 277714 + 2777i 5 + 2777i 6 + p 1 - £> 4 

-P5 -P8 -P9 -Pl2 -Pl4 -P15), (31) 

777 5 = ^(1 + 2777 2 + 2777 10 + 277714 + 2t77 15 + 2777i 6 - p 1 - p A 

+ P5-P8-P9-P12-PU-P15), (32) 

777 6 = -777 2 - 777 10 - 777i 4 + p 14 , (33) 

777 8 = -777 7 - 777i5 ~ ™16 + Pl2, (34) 

777 9 = ^(1 + 2777 3 + 2777 7 + 277714 + 277715 + 2777i 6 - p 1 - p A 

-P5-P8+P9-P12-PU-P15), (35) 

777n = -777 3 - 777 7 - 777i5 + P15, (36) 
m 12 = — mio - 77714 - "7i 6 + ps, (37) 
77713 = -m U - m 15 - 777 16 + p 4 . (38) 



By inserting the 777j's given by Eqs. (30)-(38) into Eqs. (1)-(16), and recalling the 
relations in Eq. (27), one can reproduce whichever consistent set of probabilities 
{PiiP2i ■ ■ ■ iPi&\- Indeed, as may easily be verified, the sum mi + 777 2 + 7773 + 7774, 
with mi and 7774 given by Eqs. (30) and (31), respectively, gives p\. Likewise, the 
sum 7775 + 777 6 + 7777 + 777 8 , with 7775, 777 6 , and 777 8 given by Eqs. (32), (33), and (34), 
respectively, gives p 2 = (1 -pi -p 4 +P5 ~Ps ~P% +P12 +Pu -Pis)/2, etc. Similarly, 
from Eqs. (31)-(38), we may also check that the sum 7774 + 7775 + 777 6 + 777 8 + 7779 + 
777n + m i2 + 777i3 amounts to 

Si = ^(3 - pi - p 4 - p 5 - Ps - P9 ~ P12 - Pu - Pis)- (39) 

Hence, from Eqs. (24) and (39), we obtain the LHV prediction A LH v = 2(pi +Pa + 
Ph + Ps + P9 + P12 + Pu + P15 — 2), thereby reproducing the generic Bell-CHSH sum 
of correlations in Eq. (23). Of course, according to our previous result following 
Eq. (24), a violation of the Bell-CHSH inequality necessarily implies the quantity 
Si in Eq. (39) to be either Si < or Ei > 1. We note the remarkable fact that 
there remain seven degrees of freedom in solutions for LHV probability measures 
reproducing any consistent set of probabilities. For the general solution displayed 
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in Eqs. (30)-(38), these degrees of freedom correspond to the variables rri2, m 3 , m 7 , 
m io, m u, mis, and m 16 . 

We conclude by discussing the case where, due to perfect correlation between 
the particles, two of the probabilities, say p 2 and p 3 , are equal to zero. This means 
that the results for the joint measurement of the observables a\ and b\ must both 
be either +1 or —1. Thus, from a physical point of view, it is reasonable to suppose 
that, for the case in which p 2 = and p 3 = 0, the probability measures ms, m 6 , m 7 , 
m 8 , m 9 , m 10 , m n , and m 12 do equally vanish (see Eqs. (2) and (3)). Otherwise, the 
LHV model could yield joint detection events which, by assumption, never happen. 
(Of course, mathematically, we may have p 2 = without actually requiring that 
m 5 — m 6 — m 7 — m 8 — 0- This will happen, for example, whenever m 5 + itiq = 
— m 7 — m 8 .) On the other hand, the fact that p 2 = p 3 = imposes further constraints 
on the probabilities. Specifically, since p\ +p 2 +p 3 + p 4 = 1, we have that pi = l—p4 
whenever p 2 = p 3 = 0. In addition to this, the first equation in (27) tells us that 
P5 = Ps + P9 - P12 - Pu + P15 whenever p x — 1 - p 4 and p 2 = 0. Then, putting 
1^5 = 1^6 — m 7 — m s — m 9 — m io — m n — m i2 = 0, and substituting 1 — p^ for 
Pi, and ps +P9 — Pn — Pu + P15 for ps in the general equations (30)- (38), we obtain 
the following solution for the remaining LHV probability measures in the case that 
perfect correlation develops between the measurement outcomes of the observables 
ai and bi, 



mi 


= -mi 6 +p 8 +p 9 -pu, 


m 2 


= mi 6 - p 8 +J914, 


m 3 


= mi 6 -pi2 +P15, 


m 4 


= 1 - mis - p 4 - p 9 + p 12 - pi 5 , 


m 13 


= mi 6 +p 4 -Ps ~Pl2, 


m i4 


= -m 16 +p 8 , 


mi 5 


= -mi 6 +P12, 



where now there remains a degree of freedom in the solution corresponding to the 
variable m 16 . This degree of freedom that remains in the solution was already 
noted by Han et al. [7]. The solution in Eq. (40) gives the prediction Alhv — 
2(2p§ + 2pg + 2pi5 — 1), so the Bell-CHSH inequality will be violated whenever 
Ps + P9 + P15 > 1- Incidentally, we can easily prove the negativity of either m 4 or 
mi3 in this case by simply noting that m 4 + mi 3 = 1 — p 8 — p 9 — pi 5 . 

In summary, in this paper we have proved the necessity of negative probability 
measures for the considered LHV model in all instances where the predictions by 
such a model are to give a violation of the Bell-CHSH inequality. Moreover, we 
have obtained the most general solution for LHV probability measures that repro- 
duce any conceivable set of probabilities fulfilling the normalization condition and 
the causal communication constraint. We have observed that there remain seven 
degrees of freedom in the solution. In this respect, it should be emphasized that, as 
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we have seen, it is only by imposing the condition of perfect correlation and fixing 
the eight corresponding LHV probability measures to zero, that the number of de- 
grees of freedom remaining in the solution reduces to one. In general, however, the 
solution contains more than one degree of freedom. We remark that the achieved 
general solution for LHV probability measures can be used, specifically, to repro- 
duce whichever quantum mecanical predictions for the probabilities pi,p 2 , . . . ,pw 
since, as was mentioned previously, the quantum theoretic predictions satisfy the 
consistency conditions of normalization and causal communication (see Ref. [13] and 
references therein). 

Finally we note that, by taking a somewhat different approach, Rothman and 
Sudarshan [15] (see also Ref. [16]) have arrived at essentially the same results as those 
reached in this paper. Specifically, they show that the CHSH sum of correlations 
is derivable from a master probability distribution involving 16 joint probabilities 
for four "simultaneous" spin measurements along four axes (see Eq. (4.2) of [15]), 
and that the CHSH inequality is violated if the probabilities are allowed to become 
negative. They also give explicit solutions for the 16 four-probabilities, P(+ + ++), 

P(+ + H — ),..., P( ) (see Tables II and III of [15]), that reproduce the 

standard quantum mechanical predictions.^ 

As cleverly anticipated by Feynman [17,18] (see also Ref. [11]), the important 
point arising from the discussions at hand can be summarized by saying that [15], 
"The only difference between the classical and quantum cases is that in the former 
we assume the probabilities are positive-definite." Be that as it were, I would like 
to end this paper by quoting the final sentence ending the review paper in Ref. [3]: 
"Kolmogorov's axiom may hold or not; the probability for the existence of negative 
probabilities is not negative" (italics in the original). 

Acknowledgements — I am grateful to Dimiter G. Chakalov for bringing my 
attention to the paper in Ref. [15]. 



1 For ease of comparison with the work of Rothman and Sudarshan in Ref. [15], here we write 
down the translation between their notation and ours for the 16 four- probabilities: P(+ + ++) = 
mi, P(+ + +-) = m 2 , P(+ - ++) = m 3 , P(+ - +-) = m 4 , P(+ + -+) = m 5 , P(+ + — ) = m 6 , 

P(+ - -+) = m r , P(+ ) = m 8 , P(- + ++) = m 9 , P(- + +-) = mio, P(- - ++) = 

mil, P( 1 — ) = mi2, P( — I h) = mi3, P( — I ) = mu, P( h) = mi5, and 

P( ) = mi6- From Table III of [15] we can check that, for a polarizer angle of 45° (where 

9 = 45° is just the angle giving A = 2y/2), we have P(H 1 — ) = P(+ H h) = P(+ H ) = 

P(+ ) = P(- + ++) = P(- - ++) = P(- - +-) = P(- + -+) = (1 - V2)/16 and 

P(+ + ++) = P(+ + +-) = P(+ - ++) = P(+ - -+) = P(- + +-) = P(- + ) = 

P( h) = P( ) = (1 + \/2)/16, in accordance with Eqs. (25) and (26). 
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